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ABSTRACT 

We present Ha maps at Ikpc spatial resolution for star-forming galaxies at z ^ 1, made possible by the 
WFC3 grism on HST. Employing this capability over all hve 3D-HST/CANDELS fields provides a sample of 
2676 galaxies enabling a division into subsamples based on stellar mass and star formation rate. By creating 
deep stacked Ha images, we reach surface brightness limits of 1 x 10“*^ergs“'cm“^arcsec“^, allowing us to 
map the distribution of ionized gas out to greater than lOkpc for typical L* galaxies at this epoch. We find 
that the spatial extent of the Ha distribution increases with stellar mass as r^a = 1.5(M*/1O"'M0)'''^^ kpc. 
Furthermore, the Ha emission is more extended than the stellar continuum emission, consistent with inside-out 
assembly of galactic disks. This effect, however, is mass dependent with thq/c* = I.^M^/IO'^Mq)""^^, such 
that at low masses rua ~ f*- We map the Ha distribution as a function of SFR(IRh-UV) and find evidence 
for ‘coherent star formation’ across the SFR-M* plane: above the main sequence, Ha is enhanced at all radii; 
below the main sequence, Ha is depressed at all radii. This suggests that at all masses the physical processes 
driving the enhancement or suppression of star formation act throughout the disks of galaxies. It also confirms 
that the scatter in the star forming main sequence is real and caused by variations in the star formation rate at 
fixed mass. At high masses (lO*"'^ < M,,/Mq < 10"), above the main sequence, Ha is particularly enhanced 
in the center, indicating that gas is being funneled to the central regions of these galaxies to build bulges and/or 
supermassive black holes. Below the main sequence, the star forming disks are more compact and a strong 
central dip in the EW(Ha), and the inferred specific star formation rate, appears. Importantly though, across 
the entirety of the SFR-M* plane we probe, the absolute star formation rate as traced by Ha is always centrally 
peaked, even in galaxies below the main sequence. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift — galaxies: structure 
— galaxies: star formation 


1. INTRODUCTION 

The structural formation history of galaxies is written by 
the spatial distribution of their star formation through cosmic 
time. Recently, the combination of empirical modeling and 
observations of the scaling relation between stellar mass and 
star formation rate has enabled us to constrain the build up of 
stellar mass in galaxies over a large fraction of cosmic time 
(Yang et al. 2012; Leja et al. 2013; Behroozi et al. 2013; 
Moster, Naab, & White 2013; Lu et al. 2014; Whitaker et al. 
2014). The dawn of Wide Field Camera 3 (WFC3) on the 
Hubble Space Telescope (HST) has enabled us to map the 
structural growth of this stellar mass content of galaxies at 
high fidelity over a large fraction of the history of the uni¬ 
verse (e.g. Wuyts et al. 2011a, 2012; van der Wei et al. 
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2012; van der Wei et al. 2014a, 2014b; Bruce et al. 2014; 
Boada et al. 2015; Peth et al. 2015). It has become clear that 
the physical sizes of galaxies increase with cosmic time as 
the universe expands (Giavalisco, Steidel, & Macchetto 1996; 
Ferguson et al. 2004; Oesch et al. 2010; Mosleh et al. 2012; 
Trujillo et al. 2006; Franx et al. 2008; Williams et al. 2010; 
Toft et al. 2007; Buitrago et al. 2008; Kriek et al. 2009; 
van der Wei et al. 2014a). For star forming galaxies, with 
increasing stellar mass, the disk scale length increases as does 
the prominence of the bulge (e.g. Shen et al. 2003; Lang 
et al. 2014). The picture that has emerged from these studies 
is that most galaxies form their stars in disks growing inside 
out (Wuyts et al. 2011a, 2013; van der Wei et al. 2014b; 
Abramson et al. 2014). 

In the canonical paradigm, inside-out growth is a conse¬ 
quence of the dark mater halo properties of the galaxies. 
Galaxies are thought to accrete their gas from the cosmic web 
at a rate throttled by the mass of their dark matter halo (e.g. 
White & Rees 1978; Dekel et al. 2013). The gas cools onto 
the disk of the galaxy and forms stars with a radial distribution 
set by the angular momentum distribution of the halo (Fall 
& Efstathiou 1980; Dalcanton, Spergel, & Summers 1997; 
van den Bosch 2001). As the scale factor of the universe 
increases, so does the spatial extent of the gas (Mo, Mao, 
& White 1998); galaxies were smaller in the past and grow 
larger with time, building up from the inside-out. However, 
the actual formation of galaxies in a cosmological context is 
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more complex (e.g., van den Bosch 2001; Hummels & Bryan 
2012). Recently, significant progress has been made by the 
creation of realistic disk galaxies in hydrodynamical simula¬ 
tions (Governato et al. 2010; Agertz, Teyssier, & Moore 2011; 
Guedes et al. 2011; Brooks et al. 2011; Stinson et al. 2013; 
Aumer et al. 2013; Marinacci, Pakmor, & Springel 2013) 
and combining theory and observations in a self-consistent 
framework (Keres et al. 2009; Dekel & Birnboim 2006; Dekel 
et al. 2009b; Genzel et al. 2008, 2011; Forster Schreiber et al. 
2009, 2011a; Wuyts et al. 201 lb, 2011a). How gas is accreted 
on to galaxies (e.g. Brooks et al. 2009; Sales et al. 2012) and 
feedback (e.g. Keres et al. 2005; Sales et al. 2010; Ubler et al. 
2014; Nelson et al. 2015; Genel et al. 2015) have been shown 
to be essential ingredients. However, precisely what physical 
processes drive the sizes, morphologies, and evolution of disk 
galaxies is still a matter of much debate (see, e.g., Dutton & 
van den Bosch 2012; Scannapieco et al. 2012). 

Furthermore, the evidence for this picture is indirect: we 
do not actually observe star formation building up different 
parts of these galaxies. Instead, we infer it based on empiri¬ 
cally linking galaxies across cosmic time and tracking radial 
changes in stellar surface densities and structural parameters 
(van Dokkum et al. 2010; Wuyts et al. 2011a; van Dokkum 
et al. 2013; Patel et al. 2013; van der Wei et al. 2014a; Bren¬ 
nan et al. 2015; Papovich et al. 2015). However, this method 
has uncertainties due to scatter in stellar mass growth rates 
and merging (e.g. Leja et al. 2013; Behroozi et al. 2013). Fur¬ 
thermore, migration and secular evolution may have changed 
the orbits of stars after their formation such that they no longer 
live in their birthplaces (e.g., Roskar et al. 2008). 

The missing piece is a direct measurement of the spatial 
distribution of star formation within galaxies. This is crucial 
to understanding the integrated relations of galaxy growth be¬ 
tween SFR and M*. The spatial distribution of star formation 
yields insights into what processes drive the star formation 
activity, evolution of stellar mass, and the relation between 
them. It helps to disentangle the role of gas accretion, merg¬ 
ers, and secular evolution on the assembly history of galaxies. 
Furthermore, this provides a test of inside-out growth which 
appears to be a crucial feature of galaxy assembly history. 

What is required is high spatial resolution maps of star for¬ 
mation and stellar continuum emission for large samples of 
galaxies while they were actively forming their disks. The 
Ha flux scales with the quantity of ionizing photons pro¬ 
duced by hot young stars, serving as an excellent probe of the 
sites of ongoing star formation activity (Kennicutt 1998). A 
number of large surveys have used Ha to probe the growth of 
evolving galaxies, including recently: HiZELS (Geach et al. 
2008; Sobral et al. 2009), WISP (Atek et al. 2010), MASSIV 
(Contini et al. 2012), SINS/zC-SINF (Forster Schreiber et al. 
2006, 2009), KROSS, Stott et al. (2014), and KMOS^" (Wis- 
nioski et al. 2015). Broadband rest-frame optical imaging 
provides information on the stellar component. The spatial 
distribution of this stellar light contains a record of past dy¬ 
namical processes and the history of star formation. The com¬ 
parison of the spatial distribution of ionized gas and stellar 
continuum emission thus provides an essential lever arm for 
constraining the structural assembly of galaxies. This potent 
combination shed light on the turbulent early phase of mas¬ 
sive galaxy growth at z ~ 2 (Forster Schreiber et al. 2011a; 
Genzel et al. 2014a; Tacchella et al. 2015b, 2015a), and the 
spatially-resolved star-forming sequence (Wuyts et al. 2013). 
To apply this same methodology to a global structural analysis 


requires high spatial resolution spectroscopic measurements 
for a large sample of galaxies. An ideal dataset would also 
contain broadband optical imaging with the same high spatial 
resolution to allow for robust comparison of the spatial distri¬ 
bution of ionized gas and stellar continuum emission. 

This has now become possible with the WFC3 grism ca¬ 
pability on HST. The combination of WFC3’s high spatial 
resolution and the grism’s low spectral resolution provides 
spatially resolved spectroscopy. Because this spectrograph is 
slitless, it provides a spectrum for every object in its field of 
view. This means that for every object its field of view and 
wavelength coverage, the grism can be used to create a high 
spatial resolution emission line map. The 3D-HST legacy pro¬ 
gram utilizes this powerful feature for a 248 orbit NIR imag¬ 
ing and grism spectroscopic survey over the five CANDELS 
fields (van Dokkum et al. 2011; Brammer et al. 2012a, Mom- 
cheva et al. in prep). In this paper, we use data from the 
3D-HST survey to map the spatial distribution of Ha emis¬ 
sion (a tracer of star formation) and //fhoiv stellar continuum 
emission (rest-frame 7000A, a proxy for the stellar mass) for 
a sample of 2676 galaxies at 0.7<z<1.5. The Ha and stellar 
continuum are resolved on scales of 0.13". This represents 
the largest survey to date of the spatially resolved properties 
of the Ha distribution in galaxies at any epoch. This spatial 
resolution, corresponding to ~ 1 kpc, is necessary for struc¬ 
tural analysis and only possible from the ground with adaptive 
optics assisted observations on 10m class telescopes. This 
dataset hence provides a link between the high spatial res¬ 
olution imaging datasets of large samples of galaxies with 
HST and high spatial resolution emission line maps of nec¬ 
essarily small samples with AO on large ground-based tele¬ 
scopes. This study complements the large MOSDEE (Kriek 
et al. 2015) and KMOS^^ (Wisnioski et al. 2015) spectro¬ 
scopic surveys by providing higher spatial resolution emission 
line measurements. 

We present the average surface brightness profiles of Ha 
and stellar continuum emission in galaxies during the epoch 
0.7 < z < 1.5. We analyze Ha maps for 2676 galaxies from 
the 3D-HST survey to trace the spatial distribution of star for¬ 
mation. Our sample cuts a large swath through the SER-M* 
plane covering two orders of magnitude in stellar mass 10^ < 
M* < lO" and star formation rate 1 < SFR < 4OOM0/yr and 
encompassing the star forming “main sequence” (MS). Wuyts 
et al. (2012) showed that the bright, visually striking clumps 
of star formation which appear to be common in high red- 
shift galaxies are short-lived and contribute little to the inte¬ 
grated SER of a galaxy. Here, we average over these short¬ 
lived clumps by stacking Ha maps. Stacking thousands of 
HST orbits provides deep average Ha images that allow us to 
trace the Ha distribution down to a surface brightness limit of 
1 X 10“'*ergs“' cm“^arcsec“^ in our deepest stacks, an order 
of magnitude fainter than previous studies in the high redshift 
universe. This enables us to measure the star formation sur¬ 
face density down to a limit of is 4 x lO-'^Meyr-'kpc-^ With 
these deep stacked images, the primary goals of this study are 
to derive the average surface brightness profile and effective 
radius of Ha as a function of mass and star formation rate to 
provide insight into where star formation occurs in galaxies at 
this epoch. 
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Fig. 1.— Sample selection. The left panel shows the location of our flux cut with respect to the locus of galaxies in the F(Ha)-M* plane. The right panel shows 
the distribution of our sample in the SFR-M* plane. The SFRs come from the UV+IR. The parent sample is shown in gray; selected galaxies are shown in black. 
The fraction of the total parent sample above the Fla flux limit and extraction magnitude limit are listed at the bottom in gray. As expected, we are significantly 
less complete at low masses and star formation rates. Further, about one third of galaxies are removed due to contamination of their spectra by other sources in 
the field. Of the galaxies above the flux and extraction limits, the fraction remaining as part of the final selection are listed in black. Our sample contains 2676 
galaxies from 0.7<z<1.5 spanning two decades in M* and SFR. 


2. DATA 

2.1. The 3D-HST Survey 

We investigate the spatial distribution of star formation in 
galaxies during the epoch spanning 0.7 < z < 1.5 across the 
SFR—M^ plane using data from the 3D-HST survey. 3D-HST 
is a 248 orbit extragalactic treasury program with HST fur¬ 
nishing NIR imaging and grism spectroscopy across a wide 
field (van Dokkum et al. 2011; Brammer et al. 2012a, Mom- 
cheva et al. in prep). HST’s G141 grism on Wide Field Cam¬ 
era 3 (WFC3) provides spatially resolved spectra of all objects 
in the field of view. The G141 grism has a wavelength range 
of 1.15/im < A < 1.65/im, covering the Ha emission line 
for 0.7 < z < 1.5. Combined with the accompanying//fi 4 ow 
imaging, 3D-HST enables us to derive the spatial distribution 
of Ha and rest-frame R-band emission with matching 1 kpc 
resolution for an objectively selected sample of galaxies. 

The program covers the well-studied CANDELS fields 
(Grogin et al. 2011; Koekemoer et al. 2011) AEGIS, COS¬ 
MOS, GOODS-S, UDS, and also includes GOODS-N (GO- 
11600, PI; B. Weiner.) The optical and NIR imaging from 
CANDELS in conjunction with the bountiful public photo¬ 
metric data from 0.3 - 24/rm provide stringent constraints on 
the spectral energy distributions (SEDs) of galaxies in these 
fields (Skelton et al. 2014). 

2.2. Determining z, M*, SFR 

This study depends on robustly determining galaxy inte¬ 
grated properties, specifically M* and SER. Both of these 
quantities in turn depend on a robust determination of red- 
shift and constraints on the spectral energy distributions of 
galaxies across the electro-magnetic spectrum. To do this, 
the photometric data was shepherded and aperture photom¬ 
etry was performed to construct psf-matched, deblended. 


Jf\ 25 w/H piAm/Hpimw selected photometric catalogs (see 
Skelton et al. 2014). These photometric catalogs form the 
scaffolding of this project upon which all the remaining data 
products rest. Eor this study, we rely on the rest-frame colors, 
stellar masses, and star formation rates. All of these quantities 
were derived based on constraints from across the electromag¬ 
netic spectrum. 

Our redshift fitting method also utilizes the photometry. 
This is probably not strictly necessary for the sample of Ha 
line emitting galaxies used for this study, although it helps 
to confirm the redshift of galaxies with only one emission 
line detected. It is crucial, however, for galaxies without sig¬ 
nificant emission or absorption features falling in the grism 
spectrum. To measure redshifts, the photometry and the two- 
dimensional G141 spectrum were fit simultaneously with a 
modified version of the EAzY code (Brammer, van Dokkum, 
& Coppi 2008). After finding the best redshift, emission line 
strengths were measured for all lines that fall in the grism 
wavelength range (see Momcheva et al. in prep). 

Galaxy stellar masses were derived using stellar population 
synthesis modeling of the photometry with the EAST code 
(Kriek et al. 2009). We used the Bruzual & Chariot (2003) 
templates with solar metallicity and a Chabrier (2003) initial 
mass function. We assumed exponentially declining star for¬ 
mation histories and the Calzetti et al. (2000) dust attenuation 
law (see Skelton et al. 2014). Errors in the stellar mass due 
to contamination of the broadband flux by emission lines are 
not expected to be significant for this study (see appendix in 
Whitaker et al. 2014). 

Galaxy star formation rates in this work were computed by 
summing unobscured (UV) plus dust absorbed and re-emitted 
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emission (IR) from young stars: 

SFR = SFRyv+zsCMgyr-') = 1.09 x \Q-^\Ljr + 2.2Luv)/Lq 

( 1 ) 

(Bell et al. 2005). Luy is the total UV luminosity from 1216 
- 3000 A. It is derived by scaling the rest-frame 2800 A lumi¬ 
nosity determined from the best-fit SED with EAzY (Bram- 
mer et al. 2008). L/r is the total IR luminosity from 8- 
1000/rm. It is derived by scaling the MIPS 24/rm flux density 
using a luminosity-independent template that is the log aver¬ 
age of the Dale & Helou (2002) templates with 1 < a < 2.5 
(Wuyts et al. 2008; Franx et al. 2008; Muzzin et al. 2010). 
See Whitaker et al. (2014) for more details. 

2.3. Sample Selection 

We consider all galaxies 1) in the redshift range 0.7 < 
z < 1.5 for which the Fla emission line falls in the G141 
grism wavelength coverage; 2) that have stellar masses 9.0 < 
log(M*) < 11.0, a mass range over which our //-band se¬ 
lected catalogs are complete; and 3) that are characterized 
as star-forming according to the UVJ-color criterion based 
on SED shape (Labbe et al. 2005; Wuyts et al. 2007; 
Whitaker et al. 2011). The UVJ selection separates quies¬ 
cent galaxies from star forming galaxies using the strength of 
the Balmer/4000A break which is sampled by the rest-frame 
U-V and V-J colors. These three criteria result in a parent 
sample of 8068 star-forming galaxies. The grism spectra are 
fit down to HfiArOw = 24, trimming the sample to 6612. 

We select galaxies based on a quite generous cut in Ha 
Elux: F(Ha) >3x10 *^erg/s/cm^. This limit corresponds to 
a median signal to noise S/N(Hq) = 2 and sample of 4314 
galaxies. Galaxies with lower Ha fluxes were removed as 
they may have larger redshift errors. We note here that this 
sample is Ha-limited, not Ha-selected. That is, it is a mass- 
selected sample of star-forming galaxies where we require an 
Ha flux to ensure only galaxies with correct redshifts are in¬ 
cluded. As a result of the flux and grism extraction limits, we 
are less complete at low masses and star formation rates. We 
exclude 178 galaxies which were flagged as having bad GAL- 
EIT (Peng et al. 2002) fits in the van der Wei et al. (2014a) 
catalogs, often indicative of oddities in the photometry. We 
identify galaxies that are likely to host active galactic nuclei 
(AGN) as sources with X-ray luminosity > 10"^^ ^erg s“' or 
Ha emission line widths of cr > 1000km/s (see next section). 
We remove these 57 galaxies from the sample as emission 
from AGN would complicate the interpretation of the mea¬ 
sured Ha distributions. 

Einally, of this sample, we discard 34% of galaxies due to 
contamination of their spectra by the spectra of other nearby 
objects (see next section for more detail). The contaminat¬ 
ing spectra are primarily bright stars and galaxies unrelated 
to the object, but it is possible that this criterion might lead 
to a slight bias against denser environments. The fraction of 
galaxies removed from the sample due to contamination does 
not vary with stellar mass or star formation rate. The final 
sample contains 2676 galaxies and is shown in Fig.[T] 


3. ANALYSIS 

3.1. Morphological Infonnation in the Spectrum 



Fig. 2.— Illustration of the creation of Hck emission line maps from HST 
WFC3 grism data. The top panel shows the 2D, interlaced grism spectrum. 
The second panel shows a model for the “contamination”: the spectra of all 
objects in the field except the object of interest. The third panel is a 2D model 
for the continuum emission of the galaxy. The bottom panel is the original 
spectrum with the contaminating emission from other obejcts, and the stellar 
continuum, subtracted. The result is a 2D map of the line emission at the 
spatial resolution of HST (see Sect. 3.2 for details). 

The Ha maps at the heart of this analysis are created from 
the two-dimensional 3D-HST grism spectra. The creation of 
Ha emission line maps is possible as a consequence of a 
unique interaction of features: WEC3 has high spatial resolu¬ 
tion (0."14) and the G141 grism has low (R~130) point source 
spectral resolution. A G141 grism spectrum is a series of high 
resolution images of a galaxy taken at 46A increments and 
placed next to each other on the WEC3 detector. An emission 
line in such a set up effectively emerges as an image of the 
galaxy in that line superimposed on the continuum. A reso¬ 
lution element for a galaxy at z ~ 1 corresponds to a velocity 
dispersion of cr ^ 1000 km/s, so a spectrum will only yield 
velocity information about a galaxy if the velocity difference 
across that galaxy is more than 1000 km/s . Few galaxies have 
such large line widths. Thus in general, structure in an emis¬ 
sion line is due to morphology, not kinematics. While in a 
typical ground based spectroscopy, the shape of the emission 
line yields spectral information, in our spectra it yields spatial 
information. The upshot of this property is that by subtract¬ 
ing the continuum from a spectrum, we obtain an emission 
line map of that galaxy. A sample G141 spectrum is shown in 
Fig.|2]and sample Ha maps are shown in Fig.[3] 

We note that although it is generally true that the spectral 
axes of these Ha maps do not contain kinematic information, 
there is one interesting exception: broad line AGN. With line 
widths of > lOOOkm/s, the spectra of these objects do contain 
kinematic information. These sources are very easy to pick 
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Fig. 3.— High resolution Ha maps for z ~ 1 galaxies from HST and their corresponding rest-frame optical images. The Ha generally follows the optical light 
but not always (see also Wuyts et al. 2013). 


out: they appear as point sources in the spatial direction and 
extended in the spectral direction. 

Furthermore, because the WFC3 camera has no slits, we get 
a 2D spectrum of every object in the camera’s field of view. 
For all galaxies with 0.7<z<1.5, that have an Ha emission 
line in G141’s wavelength coverage, we obtain an Ha map to 
the surface brightness limits. Based on our selection criteria, 
using this methodology, we have a sample of 2676 galaxies at 
0.7<z<1.5 with spatially resolved Ha information. 

3.2. Making Ha maps 

The reduction of the 3D-HST spectroscopy with the G141 
grism and imaging with the Hfi 4 m filter was done using 
a custom pipeline. HST data is typically reduced by driz¬ 
zling, but the observing strategy of 3D-HST allows images to 
be interlaced instead. With this dither pattern, four images 
are taken with pointing offsets that are multiples of half pix¬ 
els. The pixels from these four uncorrected frames are then 
placed on an output grid with 0.06" pixels (van Dokkum et al. 
2000). Interlacing improves the preservation of spatial in¬ 
formation, effectively improving the spatial resolution of the 
images. Crucially, interlacing also eliminates the correlated 
noise caused by drizzling. This correlated noise is problem¬ 
atic for analysis of spectroscopic data because it can masquer¬ 
ade as spectral features. 

Although the background levels in NIR images taken from 
space are lower than in those taken from earth, they are still 
significant. The modeling of the background in the grism data 
is complicated because it is composed of many faint higher 
order spectra. It is done using a linear combination of three 
physical eigen-backgrounds: zodiacal light, metastable He 
emission (Brammer et al. 2014), and scattered light from the 
Earth limb (Brammer et al. in prep). Residual background 
structure in the wavelength direction of the frames is fit and 
subtracted along the image columns. (For more information 
see Brammer et al. 2012a, 2014, Momcheva et al. in prep) 
The 2D spectra are extracted from the interlaced G141 frames 
around a spectral trace based on a geometrical mapping from 


the location of their F140W direct image positions. A sample 
2D spectrum and a pictorial depiction of the remainder of this 
subsection is shown in Fig. |2l 

The advantage of slitless spectroscopy is also its greatest 
challenge: flux from neighboring objects with overlapping 
traces can contaminate the spectrum of an object with flux 
that does not belong to it. We forward-model contamination 
with a flat spectrum based on the direct image positions and 
morphologies of contaminating objects. A second iteration is 
done to improve the models of bright {H < 22) sources us¬ 
ing their extracted spectra. An example of this contamination 
model is shown in the second panel of Fig.|2] (See Brammer 
et al. 2012a, 2012b, 2013, Momcheva et al. in prep). To re¬ 
move contamination from the spectra, we subtract these mod¬ 
els for all galaxies in the vicinity of the object of interest. 
Furthermore, for the present analysis, all regions predicted to 
have contamination which is greater than a third of the aver¬ 
age G141 background value were masked. This aggressive 
masking strategy was used to reduce the uncertainty in the in¬ 
terpretation of the Ha maps at large radii where uncertainties 
in the contamination model could introduce systematics. 

The continuum of a galaxy is modeled by convolving 
the best fit SED without emission lines with its combined 
Jf\ 25 w /Hfum /Hpieow image. The continuum model for our 
example galaxy is shown in the third panel of Eig.|2] This 
continuum model is subtracted from the 2D grism spectrum, 
removing the continuum emission and simultaneously cor¬ 
recting the emission line maps for stellar absorption. What 
remains for galaxies with 0.7 < z < 1.5 is a map of their 
Ha emission. Eive sample Ha maps and their correspond¬ 
ing Hfi 4 ow images are shown in Eig. [3 Crucially, the Ha 
and stellar continuum images were taken with the same cam¬ 
era under the same conditions. This means that differences 
in their spatial distributions are intrinsic, not due to differ¬ 
ences in the PSE. The spatial resolution is ^1 kpc for both the 
Hpuow stellar continuum and Ha emission line maps. 

The final postage stamps we use in this analysis are 80 x 80 
pixels. An HST pixel is 0.06", so this corresponds to 4.8 x 
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Fig. 4.— Illustration of the creation of Ho: image stack and the derivation of radial profiles. The panels on the left show four of the 377 Ho; maps that are 
summed to create the stack on the right. The stack is masked with the “double pacman" mask shown, in order to mitigate the effects of redshift uncertainties 
and [S II] AA6716,6731A (see § 3.3). The surface brightness profiles derived from this stack are shown above it. The raw profile is shown in black. The profile 
con'ected for residual continuum is shown in green and the profile corrected for the effects of the psf is shown in orange. 


4.8" or 38 x 38kpc at z ^ 1. Many of these postage stamps 
have a small residual positive background (smaller than the 
noise). To correct for this background, we compute the me¬ 
dian of all unmasked pixels in the 2 kpc edges of each stamp 
and subtract it. This means that we can reliably trace the sur¬ 
face brightness out to 17 kpc. Beyond this point, the surface 
brightness is definitionally zero. 

3.3. Stacking 

To measure the average spatial distribution of Ha during 
this epoch from z = 1.5 - 0.7, we create mean Ha images by 
stacking the Ha maps of individual galaxies with similar M* 
and/or SFR (See §4&5). Many studies first use Ha images 
of individual galaxies to measure the spatial distribution of 
star formation then describe average trends in this distribu¬ 
tion as a function of M* or SFR (e.g., Forster Schreiber et al. 
2006; Epinat et al. 2009; Forster Schreiber et al. 2009; Gen- 
zel et al. 2011; Nelson et al. 2012; Epinat et al. 2012; Contini 
et al. 2012; Wuyts et al. 2013; Genzel et al. 2014a). Instead, 
we first create average Ha images by stacking galaxies as a 
function M* and SER then measure the spatial distribution 
of star formation to describe trends. This stacking strategy 
leverages the strengths of our data; Ha maps taken under uni¬ 
form observing conditions for a large and objectively defined 
sample of galaxies. Erom a practical standpoint, the method¬ 
ology has the advantage that we do not need data with very 
high signal-to-noise. As a consequence, we can explore rela¬ 
tively uncharted regions of parameter space. In particular, we 
can measure the radial distribution of star formation in galax¬ 
ies across a vast expanse of the SFR-M* plane down to low 
masses and star formation rates. Additionally, we can probe 
the distribution of ionized gas in the outer regions of galaxies 
where star formation surface densities are thought to be very 


low. 

We created the stacked images by summing normalized, 
masked images of galaxies in Hfi 4 ow and Ha. To best control 
for the various systematics described in the remainder of this 
section, for our primary analysis, we do not distort the galaxy 
images by de-projecting, rotating, or scaling them. We show 
major-axis aligned stacks in § 6 and de-projected, radially- 
normalized profiles in an appendix. Our results remain quali¬ 
tatively consistent regardless of this methodological decision. 
Eor all analyses, the images were weighted by their T/fhoiv 
flux so the stack is not dominated by a single bright object. 
The Hfi 4 Qw filter covers the full wavelength range of the 
G141 grism encompassing the Ha emission line. Normal¬ 
izing by the Hfi 4 ow emission hence accounts for very bright 
Ha line emission without inverse signal-to-noise weighting as 
normalizing by the Ha emission would. 

As a consequence of the grism’s low spectral resolution, 
we have to account for the blending of emission lines. With 
a EWHM spectral resolution of ~ 100 A, Ha A6563A and 
[N II] A6548-I-6583A are blended. To account for the con¬ 
tamination of Ha by [N ll], we scale the measured flux down 
by a factor of 1.15 (Sanders et al. 2015) and adopt this 
quantity as the Ha flux. This is a simplistic correction as 
[N Il]/Ha varies between galaxies (e.g. Savaglio et al. 2005; 
Erb et al. 2006b; Maiolino et al. 2008; Zahid et al. 2013; 
Leja et al. 2013; Wuyts et al. 2014; Sanders et al. 2015; 
Shapley et al. 2015) as well as radially within galaxies (e.g. 
Yuan et al. 2011; Queyrel et al. 2012; Swinbank et al. 2012; 
Jones et al. 2013, 2015; Eorster Schreiber et al. 2014; Genzel 
et al. 2014b; Stott et al. 2014). Stott et al. (2014) find a range 
of metallicity gradients -0.063 < AZ/Ar < 0.073 dexkpc”*, 
with the median of ^ 0 (no gradient) for 20 typical star- 
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forming galaxies at z ~ 1. Hence, we choose to adopt a single 
correction factor so as not to introduce systematic uncertain¬ 
ties into the data. 

Additionally, HaA6563A and [S ll] AA6716,6731A are 
resolved but are separated by only ^ 3 resolution elements. In 
this study, we are concerned primarily with the radial distri¬ 
bution of Ha emission. In order to prevent [S ll] from adding 
flux at large radii, we mask the region of the 2D spectrum red- 
ward of Ha where [S ll] emission could contaminate the Ha 
maps. 

Galaxies are centered according to the light-weighted cen¬ 
ter of their //fi 40 iv flux distribution. Given that the HpxjKm 
can be used as a proxy for stellar mass, we chose to center the 
galaxies according to their Hf\ 4 Qw center as our best approxi¬ 
mation of centering them according to stellar mass. While the 
HpxAm centroid will not always be the exact center of mass, 
it is a better estimate than our other option, the Ha centroid. 
We measure the centroid of the Hpij^m images as the flux- 
weighted mean pixel in the x- and y- directions independently 
with an algorithm similar to the iraf task imcntr. We shift the 
Hp\i\()w image with sub-pixel shifts using damped sine inter¬ 
polation. The G141 image is shifted with the same shifts. To 
center the Ha map requires only a geometric mapping in the 
spatial direction of 2D grism spectrum. In the spectral direc¬ 
tion, however, the redshift of a galaxy and the spatial distri¬ 
bution of its Ha are degenerate. As a result, the uncertainty 
in the spectral direction of the Ha maps is ~ 0.5 pixels (see 
Brammer et al. 2012a). 

To simultaneously address these problems, we apply an 
asymmetric double pacman mask to the Ha maps. This mask 
is shown applied to the stack in Fig.|4] The mask serves three 
purposes. First, it masks the [S ll] emission which otherwise 
could masquerade as Ha flux at large radii. Second, it miti¬ 
gates the effect of the redshift-morphology degeneracy by re¬ 
moving the parts of the Ha distribution that would be most 
affected. Third, it reduces the impact of imperfect stellar con¬ 
tinuum subtraction by masking the portion of the spectrum 
that would be most afflicted. 

A mask was also created for each galaxy’s HpiArOw image 
to cover pixels that are potentially affected by neighboring 
objects. This mask was constructed from the 3D-HST photo¬ 
metric data products. SExtractor was run on the combined 
Jf\ 25 w/H puQw/Hpiaxw detection image (see Skelton et al. 
2014). Using the SExtractor segmentation map, we flagged 
all pixels in a postage stamp belonging to other objects and 
masked them. Eor both Ha and Hp\/\Qw a bad pixel mask is 
created for known bad or missing pixels as determined from 
the data quality extensions of the fits files. 

The final mask for each Ha image is comprised of the union 
of three separate masks: 1) the bad pixel mask, 2) the asym¬ 
metric double pacman mask, and 3) the contamination mask 
(see previous section). A final Hpuow mask is made from the 
combination of two separate masks 1) the bad pixel mask and 
2) the neighbor mask. The Ha and //fhovt images are multi¬ 
plied by these masks before they are summed. Summing the 
masks creates what is effectively a weight map for the stacks. 
The raw stacks are divided by this weight map to create the 
final exposure-corrected stacked images. 

3.4. Surface brightness profiles 

The stacked Ha image for galaxies with 10'° < M* < lO'*’ ^ 
is shown in Fig.|4] With hundreds of galaxies, this image is 
very deep and we can trace the distribution of Ha out to large 
radii (~ 10 kpc). To measure the average radial profiles of the 


Ha cindHpuow emission, we compute the surface brightness 
as a function of radius by measuring the mean flux in circular 
apertures. We checked that the total flux in the stacks matched 
the Ha and Hpij^m fluxes in our catalogs. We compute error 
bars on the radial profiles by bootstrap resampling the stacks 
and in general, we cut off the profiles when S/N <2.5. The 
Ha profile for the example stack is shown in Eig.|4] Before 
moving on to discussing the trends in the observed radial pro¬ 
files, we note two additional corrections made to them. 

First, we correct the continuum model used to create the Ha 
maps. This continuum model goes out to the edge of the seg¬ 
mentation map of each galaxy, which typically encompasses 
> 95% of the light. We subtract the remaining continuum flux 
by correcting the continuum model to have the same spatial 
distribution as the broad band light. The //fhoiv filter covers 
the same wavelength range as the G141 grism. Therefore, the 
radial distribution of T/fmovt emission reflects the true radial 
distribution of continuum emission. We derive a correction 
factor to the continuum model of each stack by fitting a sec¬ 
ond degree polynomial to the radial ratio of the Hpuow stack 
to the stacked continuum model. This continuum correction 
is < 20% at all radii in the profiles shown here. 

Second, we correct the radial profiles for the effect of the 
PSF Compared to typical ground-based observations, our 
space-based PSF is narrow and relatively stable. We model 
the PSF using Tiny Tim (Krist 1995) and interlacing the 
model PSFs in the same way as the data. The FWHM is 0.14", 
which corresponds to ^ 1 kpc at z ^ 1. Although this is small, 
it has an effect, particularly by blurring the centers of the ra¬ 
dial profiles. Images can be corrected using a deconvolution 
algorithm. However, there are complications with added noise 
in low S/N regions and no algorithm perfectly reconstructs the 
intrinsic light distribution (see e.g. van Dokkum et al. 2010). 
We instead employ the algorithmically more straight-forward 
method of Szomoru et al. (2010). This method takes advan¬ 
tage of the GALFIT code which convolves models with the 
PSF to fit galaxy light distributions (Peng et al. 2002). We 
begin by fitting the stacks with Sersic (1968) models using 
GALFIT (Peng et al. 2002). These Sersic fits are quite good 
and the images show small residuals. We use these fit param¬ 
eters to create an unconvolved model. To account for devia¬ 
tions from a perfect Sersic fit, we add the residuals to this un¬ 
convolved image. Although the residuals are still convolved 
with the PSF, this method has been shown to reconstruct the 
true flux distribution even when the galaxies are poorly fit by 
a Sersic profile (Szomoru et al. 2010). It is worth noting 
again that the residuals in these fits are small so the residual- 
correction step in this procedure is not critical to the conclu¬ 
sions of this paper. 

4. THE DISTRIBUTION OE Ha AS A EUNCTION OE 
STELLAR MASS AND RADIUS 

The structure of galaxies (e.g. Wuyts et al. 2011a; van der 
Wei et al. 2014a) and their sSERs (e.g. Whitaker et al. 2014) 
change as a function of stellar mass. This means that both 
where a galaxy is growing and how rapidly it is growing de¬ 
pend on how much stellar mass it has already assembled. In 
this section, we investigate where galaxies are building stel¬ 
lar mass by considering the average radial distribution of Ha 
emission in different mass ranges. 

To measure the average spatial distribution of Ha during 
this epoch from z = 1.5 - 0.7, we create mean Ha images by 
stacking the Ha maps of individual galaxies as described in 
§ 3.3. The stacking technique employed in this paper serves 
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Fig. 5.— The average radial distribution of Ha emission in galaxies in bins of stellar mass indicated at the top of each panel. The filled circles show the 
radial profiles measured directly from the stacked Ha images. The open circles show the profiles corrected for the effect of the PSF. The lines show the best fit 
exponentials for 0.5rs < r < 3rs to the PSF-corrected profiles. There appears to be some excess flux over a pure exponential at small and large radii. The short 
vertical lines show the corresponding Ha effective radii. 



log(M.) [MJ 

Fig. 6.— Size-mass relations for Ha {rna -M,) stellar continuum (r« — 
M*). The size of star forming disks traced by Ho increases with stellar 
mass as rua oc At low masses, rHa ~ r*, as mass increases the disk 

scale length of Ho becomes larger than the stellar continuum emission as 
rHa oc Interpreting Ho as star formation and stellar continuum 

as stellar mass, this seiwes as evidence that on average, galaxies are growing 
larger in size due to star formation. 

to increase the S/N ratio, enabling us to trace the profile of 
Ha to large radii. An obvious disadvantage is that the Ha 
distribution is known to be different for different galaxies. As 
an example, the Ha maps of the galaxies shown in Fig. |3] are 
quite diverse, displaying a range of sizes, surface densities, 
and morphologies. Additionally, star formation in the early 
universe often appears to be clumpy and stochastic. Differ¬ 
ent regions of galaxies light up with new stars for short peri¬ 
ods of time. These clumps, while visually striking, make up 
a small fraction of the total star formation at any given time. 
Only 10-15% of star formation occurs in clumps while the re¬ 
maining 85-90% of star formation occurs in a smooth disk or 
bulge component (Forster Schreiber et al. 201 lb; Wuyts et al. 


2012, 2013). Stacking Ha smoothes over the short-timescale 
stochasticity to reveal the time-averaged spatial distribution 
of star formation. 

Fig.|5]shows the radial surface brightness profiles of Ha as 
a function of stellar mass. The first and most obvious fea¬ 
ture of these profiles is that the Ha is brightest in the center 
of these galaxies: the radial surface brightness of Ha rises 
monotonically toward small radii. The average distribution of 
ionized gas is not centrally depressed or even flat, it is cen¬ 
trally peaked. This shows that there is substantial on-going 
star formation in the centers of galaxies at all masses at z 1. 

With regard to profile shape, in log(flux)-linear(radius) 
space, these profiles appear to be nearly linear indicating they 
are mostly exponential. There is a slight excess at small and 
large radii compared to an exponential profile. However, the 
profile shape is dependent on the stacking methodology: if 
the profiles are deprojected and normalized by their effective 
radius (as derived from the Hfuqw data) they are closer to 
exponential (see appendix). We do not use these normalized 
profiles as the default in the analysis, as it is difficult to ac¬ 
count for the effects of the PSF. 

We quantify the size of the ionized gas distribution in two 
ways: fitting exponential profiles and Sersic models. For sim¬ 
plicity, we measure the disk scale lengths (= rj) of the ion¬ 
ized gas by fitting the profiles with an exponential between 
0.5rj < r < 3rj. These fits are shown in Fig.|5] It is clear that 
over the region 0.5rj < r <3rs the Ha distribution is reason¬ 
ably well-approximated by an exponential. Out to Sr^, ^ 90% 
of the Ha can be accounted for by this single exponential disk 
fit. This implies that most of the Ha lies in a disk. 

The scale length of the exponential disk fits increases with 
mass from 1.3kpc for 9.0 < M* < 9.5 to 2.6kpc for 10.5 < 
M* < 11.0. With re = 1.678rj, this corresponds to effective 
(half-light) radii of 2.2 kpc and 4.4 kpc respectively. We fit the 
size-mass relation of the ionized gas disks (r/za-M*) with: 

rffa(m^)= ( 2 ) 

where m* = M*/10'°Mq. Fitting the surface bright¬ 

ness profiles in the same way shows the exponential disk scale 
lengths of the stellar continuum emission vs. the ionized gas. 
We parameterize this comparison in terms of the stellar con¬ 
tinuum size: 

r,(m*) = 1.4m°''^ 


(3) 
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Fig. 7.— Average radial surface brightness profiles of ///ri4oiv (left), Ha (center), and average radial Ha equivalent width profile (EW(Ha)) (right) in galaxies 
as a function of stellar mass. The radial EW(Ha) profile is the quotient of the Ha and stellar continuum profiles, providing a comparison between the spatial 
distribution of Ha and stellai* continuum emission. At low masses the EW(Ha) profile is flat. As mass increases EW(Ha) rises increasingly steeply from 
the center, showing, in agreement with the larger disk scale lengths of Fig. [6] that the Ha has a more extended distribution than the existing stellar continuum 
emission. 


rHa(m *, r*) = 1.1 r, (4) 

For IO^Mq < M* < Mq, the Ha emission has the same 
disk scale length as the Hf 140 ^ emission. This suggests that 
the Ha emission closely follows the Hfi 4 ow emission (or 
possibly the other way around). At stellar masses M* > 10®'^ 
the scale length of the Ha emission is larger than the Hfi 4 ow 
As mass increases, the Ha grows increasingly more extended 
and does not follow the Hf i 4 qw emission as closely. The size- 
mass relations for Ha and Hfuow are shown in Fig.|6] 

The ionized gas distributions can also be parameterized 
with Sersic profiles. We fit the observed, PSF-convolved 
stacks with Sersic models using GALFIT as described in the 
previous section. The Sersic index of each, which reflects the 
degree of curvature of the profile, is 1 < n < 2 for all mass 
bins, demonstrating that they are always disk-dominated. The 
Sersic indices and sizes measured with GALFIT are listed 
in Table 1. The sizes measured with GALFIT are similar 
to those measured using exponential disk fits and exhibit the 
same qualitative trends. 

While the bootstrap error bars for each individual method 
are very small, 2 - 4%, different methodologies result in sys¬ 
tematically different size measurements. We derive our de¬ 
fault sizes by fitting exponentials to the O.Sr^ < r < 3rs re¬ 
gion of PSF-corrected profiles. Fit the same way, sizes are 
10-20% larger when profiles are not corrected for the PSF 
Adopting slightly different fitting regions can also change the 
sizes by 10-20%. The GALFIT sizes are 3-15% larger. 
With all methods the trends described remain qualitatively the 
same. That is, the effective radius of the Ha emission is al¬ 
ways greater or equal to the effective radius of the Hfuow and 
both increase with stellar mass. 

The comparison between the radial distribution of Ha and 
Hfiww can be seen explicitly in their quotient, the radial Ha 
equivalent width (EW(Ha)) profile (Fig. |7]), indicating where 
the Ha emission is elevated and depressed relative to the 
Hfuow emission. The first and most obvious feature is that 
the normalization of equivalent width profiles decreases with 
increasing stellar mass, consistent with spatially-integrated 
results (Fumagalli et al. 2012) and the fact that sSFR de¬ 
clines with stellar mass (e.g. Whitaker et al. 2014). Addi¬ 
tionally, below a stellar mass of log(M),j < 9.5, the equiv¬ 
alent width profile is flat, at least on the scales of ~ Ikpc 


TABLE 1 

Structural Parameters 


log(M.) 


Ha 


Hfiaqw 

r.T 

iv 

n 

r.y 

re 

n 

9.0 < log(M J < 9.5 

I.O 

1.8 

1.9 

1.0 

1.8 

1.9 

9.5 <log(MJ < 10.0 

1.5 

2.7 

1.8 

1.3 

2.4 

1.9 

10.0 < log(M J < 10.5 

1.8 

3.2 

1.5 

1.6 

3.1 

1.7 

10.5 < log(M J < II.O 

2.6 

5.1 

1.7 

2.0 

3.9 

2.1 


Note. Disk scale length and effective radius in kpc and Sersic index for Ha 
and //fhoiv as a function of stellar mass. For an exponential disk (n=l), 
Te = 1.678ri. 

resolved by our data. These galaxies are growing rapidly 
across their disks. In addition to the overall normalization 
of the EW decreasing, as stellar mass increases the shape 
of the EW profile changes, its slope growing steeper. For 

9.5 < log(M,) < 10.0, EW(Ha) rises by a factor of - 1.3 
from the center to 2r(., for 10.5 < log(M*) < 11.0, it rises 
by ^ 3. At low masses, the entire disk is illuminated with 
new stars; at higher masses, the Ha is somewhat centrally 
depressed relative to the stellar continuum emission. Consis¬ 
tent with the measured size trends, the radial EW(Ha) profiles 
show that Ha has a similar distribution as the stellar contin¬ 
uum emission for 9.0 < log(M*) < 9.5; as mass increases Ha 
becomes more extended and less centrally concentrated than 
the stellar continuum emission. 

Interpreting Ha as star formation and Hfuoiv as stellar 
mass implies that star formation during the epoch 0.7 < z < 

1.5 is building galaxies from the inside-out as discussed in 
§7.3. 

5. THE RADIAL DISTRIBUTION OE Ha ACROSS THE 
STAR EORMING SEQUENCE 

In the previous section, we showed how the radial distribu¬ 
tion of star formation depends on the stellar mass of a galaxy. 
Here we show how it depends on the total star formation rate 
at fixed mass. In other words, we show how it depends on a 
galaxy’s position in the SER-M* plane with respect to the star 
forming main sequence. (The star forming ’main sequence’ is 
an observed locus of points in the SFR-M* plane Brinchmann 
et al. 2004; Zheng et al. 2007; Noeske et al. 2007; Elbaz 
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Fig. 8 .— We investigate the spatial distribution of star formation in galaxies 
across the SFR(UV+IR)-M* plane. To do this, we stack the Ho maps of 
galaxies on the star forming sequence main sequence (black) and compare 
to the spatial distribution of Ha in galaxies above (blue) and below (red) 
the main sequence. The parent sample is shown in gray. The fractions of 
the total parent sample above the Ha flux and extraction magnitude limit are 
listed at the bottom in gray. As expected, we are significantly less complete at 
low masses, below the main sequence. About one third of selected galaxies 
are thrown out of the stacks due to contamination of their spectra by other 
sources in the field. Of the galaxies above the flux and extraction limits, the 
fractions remaining as part of the the final selection are listed and shown in 
blue/black/red and respectively. 


et al. 2007; Daddi et al. 2007; Salim et al. 2007; Damen et al. 
2009; Magdis et al. 2010; Gonzalez et al. 2010; Karim et al. 
2011; Huang et al. 2012; Whitaker et al. 2012, 2014) 

5.1. Definition of the Star Forming Main Sequence 

We define the star forming sequence according to the results 
of Whitaker et al. (2014), interpolated to z = 1. The slope 
of the relation between SFR and M* decreases with M*, as 
predicted from galaxy growth rates derived from the evolution 
of the stellar mass function (Leja et al. 2015), reflecting the 
decreased efficiency of stellar mass growth at low and high 
masses. Whitaker et al. (2014) find that the observed scatter 
is a constant a = 0.34 dex with both redshift and M*. 

We investigate where ‘normal’ star-forming galaxies were 
forming their stars at this epoch by determining the radial 
distribution of Ha in galaxies on the main sequence. We 
elucidate how star formation is enhanced and suppressed in 
galaxies by determining where star formation is "added" in 
galaxies above the main sequence and "subtracted" in galax¬ 
ies below the main sequence. To determine where star for¬ 
mation is occurring in galaxies in these different regions of 
the SFR-M* plane, we stack Ha maps as a function of mass 
and SFR. We define the main sequence as galaxies with SFRs 
±1.2tT = ±0.4dex from the Whitaker et al. (2014) main se¬ 
quence line at z ~ 1. Specihcally, we consider galaxies ‘be¬ 
low’, ‘on’, or ‘above’ the star forming main sequence to be 
the regions [-0.8,-0.4] dex, [-0.4 ,h- 0.4] dex, or [h-0.4,h- 1.2] dex 
with respect to the main sequence line in the SFR-M* plane. 
To dehne these regions consistently we normalize the SFRs 
of all galaxies to z ~ 1 using the redshift evolution of the nor¬ 
malization of the star forming sequence from Whitaker et al. 


(2012). These dehnitions are shown pictorially by Fig. 0in 
red, black, and blue respectively. We imposed the H-1.2 dex 
upper limit above the main sequence so the stacks wouldn’t 
be dominated by a single, very bright galaxy. We impose 
the -0.8 dex due to the Ha flux-driven completeness limit. 
Fig. [S] also shows which galaxies were actually used in the 
stacks. Our broad band magnitude extraction limit and Ha 
flux limit manifest themselves as incompleteness primarily at 
low masses and SFRs as reflected in the gray numbers and 
hlled symbols. 

We adopted this ±\.2a definition of the main sequence to 
enable us to probe the top and bottom 10% of star formers and 
ferret out differences between galaxies growing very rapidly, 
very slowly, and those growing relatively normally. Accord¬ 
ing to our dehnition (±1.2cr), the ‘Main Sequence’ accounts 
for the vast majority of galaxy growth. It encompasses 80% 
of UVJ star-forming galaxies and 76% of star formation. The 
star forming main sequence is dehned by the running median 
star formation rate of galaxies as a function of mass. The def¬ 
inition is nearly identical when the mode is used instead, indi¬ 
cating that it defines the most common rate of growth. While 
we left 20% of star-forming galaxies to probe the extremes of 
rapid and slow growth, only 7% of these galaxies live above 
the main sequence and nearly double that, 13%, live below 
it. This is a manifestation of the fact that the distribution of 
star formation rates at a given mass is skewed toward low star 
formation rates. Counting galaxies, however, understates the 
importance of galaxies above the main sequence to galaxy 
evolution because they are building stellar mass so rapidly. 
Considering instead the contribution to the total star forma¬ 
tion budget at this epoch, galaxies above the main sequence 
account for > 20% of star formation while galaxies below the 
main sequenceonly account for < 3%. 

5.2. Results 

One of the primary results of this paper is shown in Fig.|9] 
the radial distribution of Ha on, above and below the star 
forming main sequence. Above the main sequence, Ha is el¬ 
evated at all radii. Below the main sequence, Ha is depressed 
at all radii. The prohles are remarkably similar above, on, and 
below the main sequence - a phenomenon that can be referred 
to as ‘coherent star formation’, in the sense that the offsets in 
the star formation rate are spatially-coherent. As shown in 
and Fig.[T0l the offset is roughly a factor of 2 and nearly in¬ 
dependent of radius: at r < 2kpc the mean offset is a factor 
of 2.2, at 3 < r < 5 kpc it is a factor of 2.1. Above the main 
sequence at the highest masses where we have the signal-to- 
noise to trace the Ha to large radii, we can see that the Ha 
remains enhanced by a factor of > 2 even beyond 10 kpc. The 
most robust conclusion we can draw from the radial prohles 
of Ha is that star formation from ~ 2 - 6 kpc is enhanced in 
galaxies above the main sequence and suppressed in galaxies 
below the main sequence (but see § 7.4 for further discussion). 

We emphasize that the SFRs used in this paper were derived 
from UVh-IR emission. These star formation rate indicators 
are measured independently from the Ha Hux. Thus, it is not 
a priori clear that the Ha emission is enhanced or depressed 
for galaxies above or below the star forming main sequence as 
derived from the UVh-IR emission. The fact that it is implies 
that the scatter in the star forming sequence is real and caused 
by variations in the star formation rate (see § 7.4). 

In the middle panels of Fig.|9]we show the radial prohles of 
Hfxaqw emission as a function of M* above, on, and below 
the star forming main sequence. As expected, we hnd that 
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Fig. 9.— Radial surface brightness profiles of Ha, F/fmoiv. and their ratio EW(Ha) as a function of M* and SFR. The colors delineate position with respect to 
the star forming ‘main sequence’: above (blue), on (black), and below (red). Above the star forming main sequence, the Ha (as well as the and EW(Ha) 

is elevated at all radii. Below the stai* forming main sequence, the Ha is depressed at all radii. The average radial profiles are always centrally peaked in Ha and 
never centrally peaked in EW(Ha). 


the average sizes and Sersic indices of galaxies increase with 
increasing stellar mass. Disk scale lengths of Ha and Hf 14 m 
are listed in Table 2. At high masses, we find that above and 
below the main sequence, the is somewhat more cen¬ 

trally concentrated than on the main sequence (consistent with 
Wuyts et al. 2011a; Lang et al. 2014, Whitaker et al. in prep), 
possibly indicating more dominant bulges below and above 
the main sequence. We note that these trends are less obvi¬ 
ous at lower masses. Furthermore, as one would expect, the 
mass to light ratio decreases with sSFR because young stars 
are brighter than old stars. Therefore, at fixed mass, galaxies 
above the main sequence have brighter //fuoiv stellar con¬ 
tinuum emission and galaxies below the main sequence have 
fainter//f 1401V emission. 

In the bottom panels of Fig.|9]we show the radial EW(Ha) 
profiles. The most obvious feature of these profiles is that 
EW(Ha) is never centrally peaked. EW(Ha) is always flat 
or centrally depressed, indicating the Ha is always equally 
or less centrally concentrated than the the Hf 14 m emission. 
Above the main sequence, the EW(Ha) is elevated at all radii. 
Below the main sequence, the EW(Ha) is depressed at most 
radii. These trends are discussed more extensively in § 7.4-5, 


where we convert the EW(Ha) profiles to sSER profiles. 

6. EEEECTS OE ORIENTATION 

In the previous sections we analyzed average images and 
radial profiles of Ha emission with galaxies stacked as they 
were oriented on the detector. This methodology has the ad¬ 
vantage that it allows for better control of systematics. In par¬ 
ticular, we can effectively subtract the continuum out to large 
radii as we can use the radial distribution of the Hf \ 4 qw flux 
to correct for the < 5% of flux missing from the continuum 
models. A galaxy’s position angle on the detector, however, 
is arbitrary and has no physical meaning. 

Here we present stacks of galaxies rotated to be aligned 
along the major axis, as measured from the continuum emis¬ 
sion. This is an important test of the idea that the Ha emis¬ 
sion originates in disks that are aligned with the stellar dis¬ 
tribution: in that case these rotated Ha stacks should have 
similar axis ratios as the rotated //fuoiv stacks. We divide the 
galaxies into the same mass bins as in the previous sections, 
and compare the most face-on vs. the most edge-on galax¬ 
ies. The position angle and projected axis ratio {q = B/A) of 
each galaxy is measured from its Hfuow image using GAL- 
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TABLE 2 

Disk scale lengths oe Hq and stellar continuum emission below, on, and above the star eorming main sequence 




i-jCHo) [kpc] 


rs 

(Hfuow) [kpc] 


log(M*) 

below 

MS 

above 

below 

MS 

above 

9.0 < log(M J < 9.5 

9.5 < log(M J < 10.0 
lO.O < log(MJ < 10.5 

10.5 < log(MJ <11.0 

1.43 ±0.28 

1.44 ±0.07 
1.90 ±0.14 
1 .68±0.I1 

1.24 ±0.06 
1.68 ±0.02 
1.99 ±0.05 
2.60 ±0.08 

1.12 ±0.06 
1.20±0.15 
1.95 ±0.08 

3.14 ±0.49 

1.17±0.03 

1.46 ±0.03 
1.78 ±0.08 
1.57 ±0.02 

1.24±0.01 
1.51 ±0.01 
1.83 ±0.02 
2.22 ±0.05 

1.17 ±0.03 
1.27 ±0.09 
1.82 ±0.09 
1.86±0.13 


* For an exponential disk (n=l), the half-light radius is re = 1.678ri. 



2 4 6 8 10 

r [kpc] 


Fig. 10.— Radial profiles of Ha as a function of mass normalized by the 
main sequence radial profile (MS). Above the star forming main sequence, 
the Ha is elevated at all radii (blue hues). Below the star forming main 
sequence, the Ha is depressed at all radii (red hues). 

FIT (Peng et al. 2002). We rotate the //fhow and Ha images 
according to their Hfuaw position angle to align them along 
the major axis. In each mass bin, we then create face- and 
edge-on stacks from the galaxies with the highest and lowest 
20% in projected axis ratio, respectively. 

The distribution of projected axis ratios is expected to be 
broad if most galaxies are disk-dominated (see, e.g., van der 
Wei et al. 2014b). If we interpret the galaxy images as disks 
under different orientations, we would expect the stacks of 
galaxies with the highest 20% of projected axis ratios to have 
an average axis ratio of ^ 0.9 and the stacks of galaxies with 
the lowest 20% of projected axis ratios to be flattened with 
average axis ratios of ~ 0.3 (see van der Wei et al. 2014b). 
As shown in Fig. [TT]the rotated //fhow stacks are consistent 
with this expectation. Furthermore, the rotated Ha stacks are 
qualitatively very similar to the rotated Hf\ 4 qw stacks, which 
means that the Ha emission is aligned with that of the stars. 

For the edge-on stacks, we measure the flattening of the 
Ha emission and compare it to that of the uow emission. 
In the four mass bins, from low mass to high mass, we And 
q{Ha) = [0.29 ± 0.02,0.32 ± 0.03,0.31 ± 0.02,0.37 ± 0.02] 
and ^(//fi 4 ow) = [0.28 ±0.01,0.27±0.01,0.29±0.01,0.34± 


0.01] respectively, where the errors are determined from boot¬ 
strap resampling. We And that the average axis ratio of 
Hpiww emission is = 0.295 ± 0.005 and q{Ha) = 

0.323 ±0.011. We conclude that the Ha is slightly less flat¬ 
tened than the //fmow emission, but the difference is only 
marginally significant. 

There are physical reasons why Ha can have an intrinsi¬ 
cally larger scale height than the Hfum emission. Given that 
outflows are ubiquitous in the z ^ 2 universe (e.g. Shapley 
et al. 2003; Shapiro et al. 2009; Genzel et al. 2011; New¬ 
man et al. 2012; Kornei et al. 2012; Forster Schreiber et al. 
2014; Genzel et al. 2014a), it is possible that the Ha would 
have a larger scale height due winds driving ionized gas out of 
the plane of the stellar disk. Furthermore, attenuation towards 
HII regions could be more severe in the midplane of the disk 
than outside of it. This would result in Ha emission being 
less concentrated around the plane of the disk, giving a larger 
scale height. Finally, the gas disks and the stellar disks can be 
misaligned. The fact that the edge-on Ha and //fuoiv stacks 
are so similar shows that all these effects are small. 

At a more basic level, an important implication of the simi¬ 
larity of the Ha stacks and the 7/Fi4ovr stacks is that it directly 
shows that we are not stacking noise peaks. If we were just 
stacking noise, a stack of galaxies flattened in T/fuow would 
not be flattened in Ha because the noise would not know 
about the shape of the T/fhovl emission. It is remarkable that 
this holds even for the lowest mass stack, which contains the 
galaxies with the lowest Ha S/N ratio as well as the smallest 
disk scale lengths. 

7. DISCUSSION 

Thus far, we have only discussed direct observables: 
HpiAow Ha. In this Section we explicitly interpret the 
radial profiles of Ha as radial profiles of star formation and 
the radial profiles of 7 /^hoiv as radial profiles of stellar sur¬ 
face density. 

7.1. Interpreting Ha and //fmoiv as SFR and Mass 

In §4 and §5, we showed the radial distribution of Ha, 
Hpuow^ and EW(Ha). Ha emission is typically used as a 
tracer of star formation, //f mow (rest-frame optical) emission 
as a proxy for stellar mass, and EW(Ha) for the specific star 
formation rate (sSER) (e.g. Eorster Schreiber et al. 2011a; 
Wuyts et al. 2013; Genzel et al. 2014b; Tacchella et al. 
2015b, 2015a). We do the same here to gain more physical 
insight into the observed profiles. If we assume that Ha traces 
star formation and Hfi 4 ow traces stellar mass, the profiles can 
be scaled to these physical quantities using the integrated val¬ 
ues. To derive mass surface density profiles, we ignore M/L 
gradients and apply the integrated M*/L/ri 4 ow as a constant 
scale factor at all radii. Similarly, to derive star formation sur- 
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Fig. 11.— Stacks of galaxies after rotating them so their major axes are aligned, for the 20 % of galaxies with the lowest ellipticities (“face-on”) and the 20 % 
of galaxies with the highest ellipticiticies (“edge-on”). The lowest ellipticity stacks are nearly round and the highest ellipticity stacks are highly flattened with 
a/h Si 0.3, consistent with viewing disks under different projections. The Ho stacks are remarkably similar to the 1401V stacks, demonstrating that the Ho 
emission is aligned with the Hfi 4 ow emission at all masses. 


face density profile, we ignore radial dust gradients and scale 
the Ha profiles based on the integrated SFR{UV + IR)/Lna 
ratio. The sSFR profile is then the quotient of the SFR and 
M* profiles. However, there are a number of caveats associ¬ 
ated with interpreting the //fi 4 oiv, Ha, and EW(Ha) profiles 
in this manner. 

We first assess the assumption that there are no radial gradi¬ 
ents in the SFR/Ha ratio. This assumption can be undermined 
in four ways: dust, AGN, winds, and metallicity, which have 
opposing effects. Dust will increase the SFR/Ha ratio by ob¬ 
scuring the ionizing photons from star forming regions. AGN, 
winds, and higher metallicity will reduce the SFR/Ha ratio, 
as they add ionizing photons that do not trace star formation. 


These aspects, and hence the extent to which a scaling from 
Ha to SFR is a good assumption, themselves depend on stel¬ 
lar mass and star formation rate. Dust attenuation is correlated 
with stellar mass (e.g. Reddy et al. 2006, 2010; Pannella et al. 
2009; Wuyts et al. 2011b; Whitaker et al. 2012; Momcheva 
et al. 2013). At fixed mass, dust attenuation is also correlated 
with star formation rate (Wang & Heckman 1996; Adelberger 
& Steidel 2000; Hopkins et al. 2001; Reddy et al. 2006, 
2010; Wuyts et al. 2011b; Sobral et al. 2012; Dominguez 
et al. 2013; Reddy et al. 2015). Within galaxies, dust atten¬ 
uation is anti-correlated with radius (e.g., Wuyts et al. 2012), 
as it depends on the column density. This means that SFR 
and Ha should trace each other reasonably well for low mass 
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galaxies with low star formation rates, and particularly poorly 
in the the centers of massive, rapidly star-forming galaxies 
(Nelson et al. 2014; van Dokkum et al. 2015). The same 
qualitative scalings with mass and star formation likely apply 
to the likelihood of an AGN being present, outflows, and the 
contamination of Ha by [N ll]. That is, AGN are most likely 
to haunt the centers of massive, rapidly star-forming galax¬ 
ies (e.g., Rosario et al. 2013; Forster Schreiber et al. 2014; 
Genzel et al. 2014a). [N Il]/Ha is most likely to be enhanced 
above the assumed value in the centers of massive galaxies (as 
described in §3.3). Shocks from winds may contribute to the 
Ha emission in the central regions, particularly at high masses 
(Newman et al. 2012; Forster Schreiber et al. 2014; Genzel 
et al. 2014a). The takeaway here is that we are relatively con¬ 
fident interpreting Ha as star formation at low masses, low 
SFRs, and all profiles outside of the center. We are less confi¬ 
dent for the centers of the radial profiles of massive or highly 
star-forming galaxies. 

Next, we assess the assumption that there is no radial gra¬ 
dient in the M/L ratio. Dust and AGN affect the M/L in the 
same way as SFR/Ha although less strongly (e.g. Calzetti 
et al. 2000; Wuyts et al. 2013; Marsan et al. 2015; Reddy 
et al. 2015). Galaxies growing inside-out will also have gra¬ 
dients in their stellar population ages. Since older stellar pop¬ 
ulations have higher M/L ratios, these age gradients translate 
into M/L gradients. Age and dust increase M*///f uow and 
AGN decrease it. Hence using as a proxy for M* is 

a fairly safe assumption at lower masses where age and dust 
gradients are small and AGN are rare. It is somewhat less cer¬ 
tain at high masses. We also note that the contribution of the 
Ha emission to the total ///?i 4 ow' flux is small, ^ 5%. 

As the EW(Ha) profile is the quotient of the Ha and HpxjKm 
interpreting it as a profile of sSFR is accompanied by the 
amalgam of all of the above uncertainties: dust, age, AGN, 
and metallicity. This does not necessarily mean that the sSFR 
profile is more uncertain than the profiles of star formation 
and mass, as some effects cancel. In a two component dust 
model (e.g. Calzetti, Kinney, & Storchi-Bergmann 1994; 
Chariot & Fall 2000), the light from both stars and HII re¬ 
gions is attenuated by diffuse dust in the ISM. The light from 
the HII regions is attenuated additionally by dust in the undis¬ 
sipated birth clouds. Because the continuum and line emis¬ 
sion will be affected equally by the diffuse dust, the EW(Ha) 
profile will only be affected by the extra attenuation toward 
the stellar birth clouds, not the totality of the dust column. As 
a consequence, the effect of dust on the EW(Ha) profiles is 
mitigated relative to the Ha profiles. The quantity of extra 
attenuation towards HII regions remains a matter of debate 
with estimates ranging from none (Erb et al. 2006a; Reddy 
et al. 2010) to a factor of 2.3 (Calzetti et al. 2000; Yoshikawa 
et al. 2010; Wuyts et al. 2013) and many in between (e.g., 
Eorster Schreiber et al. 2009; Wuyts et al. 2011b; Mancini 
et al. 2011; Kashino et al. 2013). As with the total attenua¬ 
tion, the quantity of extra attenuation toward HII regions ap¬ 
pears to increase with M* and SER (Price et al. 2014; Reddy 
et al. 2015). Reddy et al. (2015) And that extra attenuation 
becomes significant at SER~ 2OM0/yr. If true, extra extinc¬ 
tion should be taken into account for galaxies on the main 
sequence at the highest masses, and above the main sequence 
at Iog(M J > 9.5. The issue should be less acute for galaxies 
with low masses and SERs. The only way to definitively re¬ 
solve this question is to obtain spatially-resolved dust maps in 
the future. 


7.2. Star formation in disks 

The center panel of Eig. [12] shows the radial distribution of 
SER as a function of stellar mass derived by scaling the Ha 
profiles to the total SER(UVh-IR). The radial distribution of 
SER is consistent with being disk-dominated: as discussed 
in § 4, an exponential provides a reasonably good fit to the 
profiles and the Sersic indices are 1 < n < 2. Out to ~ 
85% of the Ha can be accounted for by a single exponential 
disk fit. Approximately 15% of the Ha emission is in excess 
above an exponential: 5% from the center (< Q.5rs) and 10% 
from large radii (> 3r^). 

Taken at face value the shape of the stacked Ha profiles 
suggests that the star formation during the epoch 0.7 < z < 1.5 
mostly happens in disks with the remainder building central 
bulges and stellar halos. In reality, of course, the universe 
is likely much more complicated. Radial dust gradients will 
make the star formation appear less centrally concentrated. 
Stacking galaxies of different sizes will make the star forma¬ 
tion appear more centrally concentrated, as shown in the ap¬ 
pendix. Additionally, the gas that Ha traces can be ionized 
by physical processes other than star formation such as AGN, 
winds, or shock heating from the halo. So with the Ha we 
observe we may also be witnessing the growth of black holes, 
excited gas being driven out of galaxies, or the shock heating 
of the inflowing gas that fuels star formation. 

Interestingly, a common feature of the Ha profiles is that 
they all peak at the center. If we interpret the Ha as star 
formation, this means that at all masses, galaxies are building 
their centers. Although we caution that shocks from winds 
and AGN could add Ha (Eorster Schreiber et al. 2014; Genzel 
et al. 2014a) and dust attenuation could subtract Ha from the 
centers of the profiles. That we observe Ha to be centrally 
peaked was not necessarily expected: recently it was found 
that some massive galaxies at z ^ 2 have Ha rings (see e.g. 
Genzel et al. 2014b; Tacchella et al. 2015a), which have been 
interpreted as evidence for inside-out quenching. We note that 
our averaged profiles do not exclude the possibility that some 
individual galaxies have rings at z ~ 1, which are offset by 
galaxies with excess emission in the center. 

7.3. Inside-out Growth 

The star formation surface density (as traced by Ha) is al¬ 
ways centrally peaked but the sSER (as traced by EW(Ha)) is 
never centrally peaked. Confirming Nelson et al. (2013) we 
And that, in general, EW(Ha) is lower in the center than at 
larger radii. Confirming Nelson et al. (2012), we And that the 
effective radius of the Ha emission is generally larger than 
the effective radius of the HpiiiQw emission. This means that 
the Ha emission is more extended and/or less centrally con¬ 
centrated than the HpiArOw emission. If Ha traces star forma¬ 
tion and Hp\/\Qw traces stellar mass, these results indicate that 
galaxies have radial gradients in their specific star formation 
rates: the sSER increases with radius. If the centers are grow¬ 
ing more slowly than the outskirts, galaxies will build out¬ 
ward, adding proportionally more stars at larger radii. This 
suggests that star formation is increasing the size of galax¬ 
ies. However, galaxies are still building significantly at their 
centers (probably even more than we see due to the effects of 
dust) consistent with the fact that size growth due to star for¬ 
mation appears to be fairly weak (van Dokkum et al. 2013; 
van der Wei et al. 2014a; van Dokkum et al. 2015). 

Additionally, there appears to be a trend in 
rs{Ha)/rs{HpuQw) with mass. Below 3 x 10 ®Mq, the 
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Fig. 12.— Stellar mass surface density(left), star formation surface density(center), and specific star formation rate (right) as a function of radius and stellar 
mass. These profiles were made by scaling the profiles of Ho, Hfnow, and EW(Hq) in Fig.|7]to SFR, M«, and sSFR. Within galaxies, the specific star formation 
rate rises radially; the star formation is more extended than the existing stellar mass. This is a direct demonstration that galaxies at this epoch grow inside-out. 


Ha and the Hfiaqw roughly trace each other: the radial EW 
profile is flat and r^iHa) ~ rs{Hf i 4 ow)- As mass increases, 
Ha becomes more extended than the Hfi^ow emission: 
the EW(Ha) profile is increasingly centrally depressed and 
rs{Ha) > rsiHfiAow)- This reflects the natural expectations 
of inside out growth and the shape of the sSER-M* relation 
from both a physical and an observational standpoint. 

Observationally, our tracers Ha and HpxjKm rnay trace 
somewhat different things as a function of increasing stellar 
mass. At the low mass end, because low mass galaxies have 
such high sSERs, it’s possible that the Hf\ 4 Qw emission is 
dominated by light from young stars and is not actually a good 
tracer of stellar mass. This means that there may in fact be a 
difference in the disk scale lengths of the stellar mass and star 
formation but it is hard to detect because our proxy for M* 
is dominated by the youngest stars. At the high mass end, 
galaxies have more dust so star formation could be preferen¬ 
tially obscured at small radii. Consequently, the Ha could ap¬ 
pear to be less centrally concentrated than the star formation 
is in reality, making the inferred size larger (see e.g. Simpson 
et al. 2015). Taken together, these effects could contribute to 
the trend of increasing rs{Ha)/rs{Hf\ 4 Qw) with stellar mass. 
However, as described in §7.1, there are a number of other 
observational effects that work in the opposite direction, de¬ 
creasing the rs{Ha) / rs{Hfi 4 ow) at high masses. Dust will also 
obscure the stellar continuum emission, meaning that the stel¬ 
lar mass could also be more concentrated than observed. Age 
gradients will also change the M/L ratio, again adding more 
stellar mass at the center. Szomoru et al. (2013) estimate 
that galaxies are ^ 25% more compact in mass than in light. 
AGN contributing line emission to the Ha profiles will also 
work to decrease this ratio by adding extra flux and decreas¬ 
ing the size of the star formation. In sum, it seems more likely 
that observational effects will increase the PsiHa)/rsiHfi 4 ow) 
with mass (and generally) than decrease it but as the effects 
act in both directions we cannot say with certainty which are 
more important. 

While many observational effects could contribute to the 
the mass dependence of the size ratio, this effect may also 
have a physical explanation. More massive galaxies have 
older mean ages. This means that a larger fraction of their 
star formation took place at earlier cosmic times. Hence, it 
is perhaps then reasonable that their stellar mass - the in¬ 
tegral of their past star formation history - would be more 


compact than the gas disks with ongoing star formation. On 
the other hand, low mass galaxies have younger mean ages, 
which means their mass-weighted sizes are closer to the sizes 
of their star forming disks. 

7.4. Above and Below the Main Sequence 

Here we return to the profiles above and below the star¬ 
forming main sequence, that is, for galaxies with relatively 
high and relatively low star formation rates for their stel¬ 
lar mass. Whitaker et al. (2012) showed that the SEDs of 
galaxies above and below the main sequence are different 
from those on it. Above the main sequence, the SEDs are 
dusty but blue which they interpreted as indicative of AGN 
or merger-induced starbursts. Below the main sequences, the 
SEDs are not dusty but red, which they interpreted as indica¬ 
tive of star formation being shut down. Additionally, Wuyts 
et al. (2011a) showed that galaxies above and below the main 
sequence were structurally more compact and centrally con¬ 
centrated than galaxies on the star forming main sequence. 

Hints as to what physical processes are driving a galaxy 
above or below the main sequence are given by these trends 
in stellar structure and SED shape. The next key piece of in¬ 
formation is where the star formation is enhanced above and 
suppressed below the star forming sequence, which we show 
here. Eor instance, if the primary physical processes driving 
galaxies above the main sequence are AGN or central star- 
bursts, we would expect Ha to be enhanced in the center but 
not at larger radii. If quenching is driven by processes acting 
from the center and progressing from the inside outward, we 
would galaxies below the main sequence to have a decrease 
in Ha primarily in the center. 

We characterize galaxies with respect to the star formation 
main sequence using their total SER(IRh-UV)s which reflect 
the total obscured-Funobscured ionizing flux from young stars. 
As described in §5.2, we find that above the main sequence, 
the Ha is enhanced at all radii; below the above the main se¬ 
quence the Ha is depressed at all radii. In Eig.[T3]we show 
SER, M*, and sSER profiles made by scaling our Ha profiles 
using the integrated SER(IRH-UV)/SER(Ha) and //fhoiv pro¬ 
files using the integrated M*/L/ri 4 oiv with all the associated 
caveats described in §7.1. Because the integrated M*/L/ri 4 oiv 
decreases with increasing SER at fixed mass, the offset in the 
Hfiaqw light profiles shown in middle panels of Eig.|9] disap¬ 
pears in the M* profiles shown middle panels of Eig.[T3] At 
fixed mass, galaxies are brighter above the main sequence and 
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Fig. 13.— Radial surface brightness profiles of SFR, M*, and their ratio sSFR as a function of M* and SFR. The colors delineate position with respect to the 
stai' forming ‘main sequence’: above (blue), on (black), and below (red). SFR and M* profiles are Ha and ///ri4ow profiles scaled to the total SFR(UV+IR) and 
M* with all the caveats described in §8.1. Out to distances as great as 8 kpc from the galactic center, star formation is enhanced in galaxies above and depressed in 
galaxies below the star forming main sequence. This is also true of the specific star formation rate. In general, the radial distribution of M* is similar on, above, 
and below the main sequence on average. It becomes slightly more centrally concentrated in galaxies above and below the main sequence at the highest masses, 
as shown in Fig. \M\ There are two take home messages from this figure: 1. The star formation rate, on average, is always the highest in the centers of galaxies. 
2. The radial distribution of star formation depends more strongly on M* than SFR at fixed mass (a galaxy’s position with respect to the main sequence). 


fainter below but the underlying mass profiles are fairly sim¬ 
ilar at all SFRs (although see next section for a discussion of 
the highest masses). On the other hand, the dust attenuation 
increases with increasing SFR at fixed mass. Acting in con¬ 
cert, dust and age mean that the EW(Ha) prohles shown in 
the bottom panels of Fig. |9] likely underestimate the true dif¬ 
ference in sSFR above, on, and below the main sequence. In 
the bottom panels of Fig.[T3the trends in sSFR are enhanced 
after accounting for dust and age. 

The most robust conclusion we can draw about the radial 
distribution of star formation, an inferred quantity, is that star 
formation in the disk between 2-6kpc is enhanced above 
the main sequence and suppressed below the main sequence. 
This, in turn, has several important implications. 

First, our results constrain the importance of AGN emis¬ 
sion above the main sequence. One possibility is that galax¬ 
ies above the star forming main sequence are there because 
the bright UV-rlR emission of an ACjN was incorrectly inter¬ 
preted as star formation. In this case, the Ha emission would 


be elevated in the center but the same as on the main sequence 
throughout the rest of the disk. This, however, is not what we 
observe; the Ha in the disk from 2-6 kpc is elevated, mean¬ 
ing that galaxies are not only above the main sequence due to 
misinterpreted AGN. 

Second, because Ha is an independent indicator of star for¬ 
mation, the fact that it is enhanced at all radii confirms that 
the scatter in the main sequence is real and due to variations 
in the star formation rate at fixed mass. If the observed main 
sequence scatter were due exclusively to measurement etTors 
in the UV-rlR SFRs, the Ha should not be enhanced or de¬ 
pressed in concert, but it is. 

Third, the prohles provide information on the importance 
of mergers and galaxy encounters “pushing” galaxies above 
the main sequence. It is well established that interaction- 
driven gravitational torques can funnel gas to the center of 
a galaxy inducing a burst of star formation (e.g., Hernquist 
1989; Barnes & Hernquist 1991, 1996; Mihos & Hernquist 
1996). However, in idealized merger simulations, Moreno 
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et al. (2015) show that while star formation is enhanced in 
the central kpc of interacting galaxies, it is suppressed ev¬ 
erywhere else. This is not what we observe: the Ha in our 
stacks above the main sequence is enhanced at all radii; it 
is not enhanced in the central kpc and suppressed at larger 
galacto-centric radii. Some ambiguity is inherent in the in¬ 
terpretation of an average distribution of Ha because the dis¬ 
tribution of Ha in individual galaxies could vary significantly 
from the average. Our stacking method cannot distinguish be¬ 
tween local enhancements at random locations in the disk and 
global enhancement of the disks of individual galaxies. Nev¬ 
ertheless, our uniformly higher star formation rates suggest 
that major mergers are not the only physical process driving 
the elevated star formation in galaxies above the star forming 
main sequence. 

Below the main sequence, it is possible the dominant pro¬ 
cesses suppressing star formation act primarily in the cen¬ 
ters of galaxies where AGN live, bulges grow, and timescales 
are short. If this were the case, we would expect Ha to be 
lower in the center of the galaxies but unchanged at large 
radii. Again, this is not what we observe: below the main se¬ 
quence Ha is suppressed at all radii, indicating that the phys¬ 
ical mechanisms suppressing star formation must act over the 
whole disk, not exclusively the center. 

Instead perhaps, for stellar masses below M - 3 X IO'^Mq, 
some cosmological hydrodynamic simulations (Sparre et al. 
2015) and models (Dutton, van den Bosch, & Dekel 2010; 
Kelson 2014) have suggested that a galaxy’s position in the 
SFR-M* plane is driven by its mass accretion history. In this 
schema, galaxies living below the main sequence had early 
formation histories and galaxies above the main sequence had 
later formation histories. Sparre et al. (2015) show that in II- 
lustris, most of the scatter in the star-forming main sequence 
is driven by these long scale {> 500Myr) features of galax¬ 
ies’ formation trajectories rather than short-term stochasticity. 
Dutton et al. (2010) predict based on this model for main se¬ 
quence scatter that the size of gas disks should be the same 
above and below the main sequence. Consistent with this pre¬ 
diction, for masses below M ~ 10'** ^M 0 , we do not see sig¬ 
nificant differences in Ha sizes above and below the main se¬ 
quence, although the error bars are large (see Table 2 for val¬ 
ues). The fact that the average radial distribution of Ha does 
not have wildly different structure above and below the main 
sequence perhaps makes more sense in the context of scat¬ 
ter driven by longer timescale variations in the mass accretion 
history as opposed to some ubiquitous physical process. In 
other words, the similarity of the radial profiles appears con¬ 
sistent with a simple model in which the overall star formation 
rate scales with the gas accretion rate (averaged over some 
timescale) and the gas distributes itself in similar structures 
regardless of its accretion rate. It will be interesting to com¬ 
pare the observed gas distributions directly to those in galaxy 
formation models. 

Regardless of the physical reasons, across the SFR-M* 
plane two important features are consistent. 1) The observed 
Ha distribution is always centrally peaked. 2) The observed 
EW(Ha) is never centrally peaked. 

7.5. Bulge growth and quenching at high masses? 

While Ha is enhanced at all radii in galaxies above the main 
sequence and suppressed at all radii below the main sequence, 
in the high mass bin (M = 10'° 10"Mq), the trends appear 

to have some radial dependence as well. To examine trends 


at high masses in more detail, in Fig. [14] we show the same 
the radial profiles of Ha, and EW(Ha) above, on and 

below the main sequence as in Eigs.l9land[T3] Here we also 
normalize by the main sequence profiles to highlight differ¬ 
ences. 

Above the main sequence, there is a central excess in Ha 
emission (left panels of Eig. [Hi. The cause of this excess 
is difficult to interpret: it could be due to an AGN or extra 
star formation in the central regions or both. As mentioned 
in § 2.3, galaxies with X-ray luminosity > 10"^^ ^erg s“' or 
a very obvious broad line component are excluded from the 
analysis in this paper. The excess central Ha emission ex¬ 
ists even when galaxies hosting obvious AGN are excluded. 
However, with a very conservative cut on broad line AGN in 
which galaxies with even marginal elongation in the spectral 
direction are excluded, the central excess in Ha disappears. 
Hence, it is possible that this central enhancement is driven 
primarily by emission from AGN. If it is due to an AGN, it 
could suggest that supermassive black holes are growing in 
this region of parameter space. If it is due to star formation, it 
could indicate that bulge construction is underway, consistent 
with the growing prominence of bulges observed during this 
epoch (Lang et al. 2014). We note that because the IR/Ha in 
this bin is so high, it is likely that the excess in central ionizing 
flux (either from star formation or an AGN) would actually be 
even larger if it were not attenuated. 

If the high SERs in galaxies above the main sequence are 
fueled by elevated gas accretion rates, the disks of these galax¬ 
ies are likely to be gas-rich. In these gas-rich environments, 
it has been suggested that gravitational torques induced by 
violent disk instability could drive gas rapidly inward by vis¬ 
cous and dynamical friction (Noguchi 1999; Dekel, Sari, & 
Ceverino 2009a; Krumholz & Burkert 2010; Bournaud et al. 
2011; Genzel et al. 2011; Eorbes, Krumholz, & Burkert 2012; 
Cacciato, Dekel, & Genel 2012; Elmegreen, Zhang, & Hunter 
2012; Dekel et al. 2013; Eorbes et al. 2014). Once in the 
center this gas could fuel the bulge and/or black hole growth 
evidenced by the excess central Ha emission. 

During the epoch 0.7 < z < 1.5 in this mass range (M = 
10 io. 5 _ loiiMg) the quenched fraction roughly doubles (from 
^ 30-60%). Since the SERs of galaxies must fall below the 
main sequence on their way to quenchdom, this region of pa¬ 
rameter space would be a good place to look for hints as to 
how galaxies quench. Relative to the main sequence, the Ha 
below the main sequence appears to be depressed in the center 
(Eig.[T4|bottom left). The Ha profile also appears depressed 
relative to the main sequence at larger radii, a manifestation 
of its smaller scale radius (Eig. [Hi. That is, we find that be¬ 
low the main sequence, the star-forming disk of Ha emission 
is both less centrally concentrated and more compact. 

In addition to the Ha in the centers of galaxies below the 
main sequence being depressed relative to galaxies on the 
main sequence, it is also depressed relative to the Hfi 4 m 
emission (Eig.[T4l top right). Interpreted as sSER, this means 
that the stellar mass doubling time in the centers of these 
galaxies is significantly lower than at larger radii. Centrally 
depressed sSER has been taken as evidence of inside-out 
quenching (Tacchella et al. 2015a). Here we show this for 
the first time explicitly below the main sequence where it is 
most straight-forward to interpret in the context of star for¬ 
mation quenching. That being said, it should be noted that 
although the Ha is centrally depressed in two interesting rel¬ 
ative senses (relative to the Hfum and relative to the main 
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Fig. 14.— Top row, radial profiles of Ha, M*, and EW(Ha) for galaxies with M = 10^^'^ —10* ^Mq. As previously, the colors delineate position with respect 
to the star forming ‘main sequence’: above (blue), on (black), and below (red). The bottom row shows these profiles normalized by the main sequence (with the 
black line divided out). Above the main sequence, the Ha (and EW(Ha)) is enhanced at all radii, but somewhat more so at small and large radii. Below the main 
sequence, the Ha (and EW(Ha)) is suppressed at all radii, but somewhat more so small and large radii. There appears to be excess central stellai* mass similaidy 
above and below the main sequence. 
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Fig. 15.— Relation between disk scale length and SFR in Ha and ///ri4ow' 
emission for galaxies with 10.5 < log(M*) < 11.0. In///rj4oiy emission, the 
disk scale length is smaller above and below the main sequence than on it. In 
Ha emission, the disk scale length below the main sequence remains smaller 
than on the main sequence but is larger above it. 


sequence Ha), in an absolute sense, the Ha is not centrally 
depressed, it is centrally peaked. That is, on average, there 
is not a hole in the observed Ha emission at the centers of 
massive galaxies below the main sequence. So while we may 
be seeing some suppression of star formation in the center of 
these galaxies below the main sequence, it is not ‘quenching’ 
in the standard sense of a complete cessation of star forma¬ 
tion. 

Our findings could be viewed in the context of an evolu¬ 
tionary pathway from bulge growth to quenching (e.g., Wuyts 
et al. 2011a; Lang et al. 2014; Genzel et al. 2014b; Tac- 
chella et al. 2015a). Consistent with Wuyts et al. (2011a), we 
hnd excess central stellar continuum emission similarly above 
and below the star forming sequence. Wuyts et al. (2011a) 
suggests that this structural similarity could indicate an evo¬ 
lutionary link between the galaxies above and below the main 
sequence. 

AGN can in principle drive gas out of the centers of their 
host galaxies, efficiently removing the fuel for star forma¬ 
tion (see e.g., Croton et al. 2006). Large bulges are also in 
principle capable of stabilizing galaxy disks and suppressing 
star formation from the inside-out (‘gravitational quenching’ 
Martig et al. 2009; Genzel et al. 2014b). Observationally, 
it seems that regardless the physical cause, galaxies quench 
after reaching a stellar surface density threshold (e.g. Franx 
et al. 2008). Whatever process is underway above the main 
sequence, there are theoretical indications that it is capable of 
suppressing star formation. Some authors argue this occurs 
from the inside-out. The deep depression in EW(Ha) in the 
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centers of galaxies below the main sequence could be taken 
as evidence for one of these quenching mechanisms acting in 
this way. 

One remaining mystery, as shown in Fig.[T5]is that the Ha 
disks have much smaller sizes below the main sequence than 
on or above it. It is possible that the galaxies below the main 
sequence formed earlier than the galaxies on or above the 
main sequence at this redshift and hence the galaxies above 
the main sequence are not actually direct progenitors of those 
below. It is also possible that these galaxies underwent some 
sort of compaction on their way to quenching (e.g. Dekel & 
Burkert 2014; Zolotov et al. 2015) The most robust thing we 
can say is that below the main sequence, Ha seems to be 
both less centrally concentrated and less extended. How ex¬ 
actly this should be interpreted is unclear without the aid of 
simulations. 

7.6. The average spatial distribution of star fonnation from 
z = 0.7-1.5 

In § 4 &5 we determined the radial profiles of star forma¬ 
tion as a function of M* and SFR. Here we briefly analyze 
the radial distribution of all star formation at this epoch, that 
is, at what distance from the center of a galaxy is a star most 
likely to form. The average Ha image of all selected galax¬ 
ies is shown in Fig. [16] This is the average spatial distribu¬ 
tion of Ha in galaxies during the epoch 0.7 < z < 1.5. Each 
galaxy has an Ha map with a depth of 2 orbits on HST. We 
summed the Ha maps of 2676 galaxies, creating the equiva¬ 
lent of a 5352 orbit Ha image. This average Ha image is 
deepest Ha image in existence for galaxies at this epoch. 
With this stacked 5352 orbit HST image, we can trace the 
radial distribution of Ha down to a surface brightness limit 
of 1 X 10 '*ergs 'em ^arcsec This allows us to map the 
distribution of Ha emission out to ~ 14kpc where the star 
formation surface density is 4 x lO-''M0yr-'kpc“2 (Kenni- 
cutt 1998). 

Weighting the radial profile of Ha by area shows its prob¬ 
ability distribution. The Ha probability distribution has a 
peak, the expectation value, at 0.75 kpc. Note, we did not 
normalize by the Hfuow flux here so the expectation value 
reflects the most likely place for a random HII region within 
a galaxy to exist. If we interpret Ha as star formation then 
during the epoch 0.7 < z < 1.5, when ^ 33% of the total star 
formation in the history of the universe occutTed, the most 
likely place for a new star to be born was 0.75 kpc from the 
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center of its home galaxy. 


8. CONCLUSIONS 

In this paper, we studied galaxy growth through star for¬ 
mation during the epoch 0.7 < z < 1.5 through a new win¬ 
dow provided by the WFC3 G141 grism on HST. This slit¬ 
less grism spectroscopy from space, with its combination of 
high spatial resolution and low spectral resolution gives spa¬ 
tially resolved Ha information, for 2676 galaxies over a large 
swath of the SFR-M* plane. Ha can be used as a proxy for 
star formation, although there are many uncertainties (§7.1). 
The most important new observational result of our study is 
the behavior of the Ha profiles above and below the main se¬ 
quence: remarkably, star formation is enhanced at all radii 
above the main sequence, and suppressed at all radii below 
the main sequence (Fig.fOl). This means that the scatter in 
the star forming sequence is real. It also suggests that the pri¬ 
mary mode of star formation is similar across all regions of 
this parameter space. 

Across the expanse of the SFR-M* plane, the radial distri¬ 
bution of star formation can be characterized in the following 
way. Most of the star formation appears to occur in disks 
(Fig.|5j, which are well-aligned with the stellar distribution 
(Fig.[IB. To first order, Ha and stellar continuum emission 
trace each other quite well. On average, the Ha surface den¬ 
sity is always highest in the centers of galaxies, just like the 
stellar mass surface density. On the other hand, the EW(Ha), 
and the inferi'ed specific star formation rate, is, on average, 
never highest in the centers of galaxies (Eig.|9j). Taken at face 
value, this means that star formation is slightly more extended 
than the existing stars (Eig.|6ll, demonstrating that galaxies at 
this epoch are growing in size due to star formation. 

The results in this study can be extended in many ways. In 
principle, the same dataset can be used to study the spatial 
distribution of [O III] emission at higher redshifts, although it 
is more difficult to interpret and the fact that it is a doublet 
poses practical difficulties. With submm interferometers such 
as NOEMA and ALMA the effects of dust obscuration can 
be mapped. Although it will be difficult to match the resolu¬ 
tion and sample size that we reach in this study, this is crucial 
as dust is the main uncertainty in the present analysis. Ei- 
nally, joint studies of the evolution of the distribution of star 
formation and the stellar mass can provide constraints on the 
importance of mergers and stellar migration in the build-up of 
present-day disks. 
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9. APPENDIX 

In this paper we investigate the average radial distribution of Ha emission by stacking the Ha maps of individual galaxies 
and computing the flux in circular apertures on this stack. With this methodology, we average over the distribution of inclination 
angles, position angles, and sizes of galaxies that go into each stack. The simplicity of this method has a number of advantages. 
First, it requires no assumptions about the intrinsic properties of galaxies. Second, it allows us to measure the average size of the 
Ha distribution in the star forming disk. Finally, because the image plane is left in tact, we can correct for the PSF. 

To complement this analysis, here we present the average deprojected, radially-normalized distribution of Ha. We do this to 
test the effect of projection and a heterogenous mix of sizes on the shape of the radial profile of Ha, to ensure trends were not 
washed out with the simpler methodology employed in the rest of the paper. 

To do this, we use GAFFIT (Peng et al. 2002) to derive the effective radius, axis ratio, and position angle of each galaxy from 
its Hf\ 4 qw stellar continuum image. We correct for the inclination angle of each galaxy by deprojecting the (x,y) pixel grid of 
it’s image based on the inclination angle implied by the axis ratio. The surface brightness profile is computed by measuring the 
flux in deprojected radial apertures. In practice, this is done simply by extracting the radial profile of each galaxy in elliptical 
apertures deflned by the position angle, axis ratio, and center of the ///?i 4 oty image. The extraction apertures were normalized 
by the Hfxaow effective radius of each galaxy. A radial profile in deprojected, r^-normalized space is derived for each galaxy. 
These individual galaxy profiles are flux-normalized by their integrated //fmow magnitude and summed to derive the mean radial 
distribution. 

The average de-projected, r^-normalized radial profiles of Ha, EW(Ha) are shown in Fig.[T2l In general, the 

qualitative trends seen here are the same as those described in the main text. For the region 0.5 < < 3 the radial profile of Ha 

remains consistent with an exponential all masses, above, on, and below the star forming sequence. The radial profiles of both Ha 
and HpiArOw are somewhat less centrally peaked than the analogous profiles in Fig.|9] This is expected of disk-dominated galaxies 
under different orientation angles as flux from the disk of edge-on galaxies could be projected onto the center. Additionally, 
stacking galaxies of different sizes can result in a somewhat steeper (higher n) profile than the individual galaxies that went into 
it (see van Dokkum et al. 2010). Because the shapes of the Ha and Hfuqw profiles are similarly effected by deriving the profiles 
with this different methodology, the shape of the EW(Ha) profiles remain largely unchanged. 
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Fig. 17.— Average de-projected, r^-normalized radial profiles of Ha, T/fmow', and EW(Ha) as a function of mass. 
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